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摘要 : 磁场 是 太阳 爆发 活动 的 根本 驱动 力 , 高 精度 实时 偏振 测量 是 实现 磁场 测量 的 常用 
方法 。 然而 常用 的 波 片 式 测量 方法 需要 进行 多 次 测量 ， 且 存在 机 械 旋 转 结构 ， 易 导致 仪器 拌 
动 ， 降 低 测量 精度 ; 波 片 旋转 或 切入 切 出 也 会 导致 测量 时 间 较 长 。 目 前 ， 基 于 DKDP 晶体 
和 波 片 的 太阳 大 气 偏振 测量 周期 较 长 , 实时 性 不 高 , 无 法 应 用 于 快速 变化 的 太阳 活动 磁场 探 
测 。 基 于 以 上 背景 ， 本 论文 提出 了 一 种 基于 双 DKDP 晶体 原理 的 太阳 大 气 实时 偏振 探测 技 
术 ， 采 用 两 片 快 轴 方 向 不 同 的 DKDP 晶体 ， 通 过 外 加 不 同 电压 进行 快速 相位 延迟 变化 ， 实 
现 对 入 射 太 阳光 的 偏振 调制 ， 获 取 太 阳 大 气 偏振 斯 托 克 斯 矢量 。 相 较 于 传统 旋转 波 片 与 


DKDP 晶体 相 结合 的 方法 ， 本 论文 提出 方法 可 将 单 次 偏振 调制 速率 从 数秒 提升 到 毫秒 级 , 极 
一 大 地 提升 了 太阳 磁场 测量 实时 性 。 仿 真 分 析 结果 表明 : 本 文 提出 方法 测量 精度 达到 2 x 10-3， 
< 单 波长 点 测量 时 间 分 辩 率 提升 了 20 倍 以 上 , 可 以 实现 对 太阳 磁场 偏振 信息 的 实时 精确 测量 
e 关键 词 ， 太 阳 磁 场 ， 偏 振 测量 ，DKDP， 实 时 测量 

= 中 图 分 类 号 : P161.3 ”文献 标识 码 ，A 
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© RRS. 这些 太阳 活动 作用 于 日 地 空间 环境 ， 对 地 球 气候 、 人 体 健康 、 电 子 电 力 设备 、 无 
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线 电 通讯 以 及 航空 飞行 器 飞行 都 有 重大 影响 中。 通过 对 太阳 活动 的 研究 来 实现 对 这 些 灾难 
性 的 空间 天 气 的 有 效 预报 , 是 当今 太阳 物理 学 研究 的 一 项 重要 课题 。 太 阳 磁 场 是 太阳 活动 的 
根本 驱动 力 ， 而 对 太阳 磁场 的 测量 通常 是 采用 对 特定 磁 敏 感 谱 线 的 侦 振 测量 来 实现 。 太 阳 谱 
线 在 磁场 中 会 发 生 分 裂 , 其 裂 距 和 文 线 侦 振 态 与 磁场 强度 和 方向 呈现 一 定 的 函数 关系 , 通过 
= 精确 测量 太阳 光 的 偏振 状态 ,建立 偏振 辐射 转移 方程 ， 反 演 太 阳 大 气 中 的 磁场 形态 ， 可 以 人 研 
FI TART REDE H eR SETS SIUS, op Has TIT. Voss TES RIAL CON 
和 预报 以 及 太阳 物理 的 研究 具有 重要 的 意义 中 。 
国内 外 对 太阳 磁场 的 观测 主要 通过 在 太阳 望远镜 上 安装 偏振 测量 仪器 对 太阳 光谱 偏振 
信息 进行 测量 实现 。 而 地 基 太 阳 望 远 镜 的 观测 会 受到 大 气 满 流 的 影响 ， 当 望远镜 放大 倍率 较 
KS, 成像 光 广 模糊 拌 动 ， 大 大 影响 观测 分 辨 力 ， 并且 白天 大 气 满 流 也 更 强 。 当 偏振 测量 的 
时 间 分 辨 率 较 高 时 , 可 在 一 定 程度 上 达到 冻结 大 气 添 流 的 效果 ,降低 大 气 漠 流 对 测量 精度 的 
影响 。 除 此 之 外 ， 太 阳 侯 场 也 始终 处 于 快速 变化 中 。 为 了 能 够 精确 测量 太阳 光 的 偏振 态 ， 反 
演 得 到 太阳 大 气 的 人 磁场， 偏振 测量 仪器 的 时 间 分 辩 率 就 成 了 一 个 重要 指标 。 
目前 , 国内 太阳 望远镜 上 使 用 的 偏振 分 析 器 主要 采用 的 是 带 有 经 典 波 片 的 调制 系统 。 这 
种 传统 的 偶 振 调制 器 也 被 许多 国外 的 太阳 侦 振 分 析 器 所 使 用 ， 如 德国 GREGOR 望远镜 上 所 
装载 的 理想 光谱 侦 振 仪 POLIS)， 其 波 片 的 旋转 频率 为 15Hz。 美 国 大 熊 湖 GST 望远镜 上 
的 可 见 光 偏振 成 像 仪 VIM》〉 和 第 二 代 近 红外 光谱 偏振 仪 “NIRIS〉 都 采用 了 旋转 波 片 调制 
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方式 。NIRIS 的 波 片 调制 器 连续 旋转 ， 其 旋转 与 红外 相机 的 帧 率 相 结合 ， 以 22.5” 步 长 旋转 
一 周 获得 16 帧 图 像 ， 目 前 其 相机 能 够 以 10 帧 / 秒 的 速度 完成 低 噪声 采集 ; 所以，NIRIS Jt 


谱 - 偏 振 测量 的 频率 通常 为 10s 
镜 一 组 8 步 式 调制 需要 4s 左 
的 时 间 与 数据 采集 时 间 之 和 ， 


;其 中 VIM 的 时 间 分 辩 率 为 10s 40, JE 
完成 ， 每 个 调制 位 置 上 耗费 的 总 时 间 为 波 片 施 转 到 位 所 需要 


间 占 了 很 大 一 部 分 54。 在 此 基础 上 ， 


而 波 片 旋转 所 消耗 的 时 


怀柔 太阳 观测 站 的 SMFT 采用 了 波 片 十 DKDP 晶体 的 调制 方式 ， 这 种 方式 可 以 利用 


内 的 NVST 望 远 


正 负电 


压 翻 转 控 制 DKDP 蝇 体 的 延迟 ， 快 速 采集 两 个 参数 的 图 像 ， 减 少 了 波 片 旋转 的 次 数 ， 缩 短 


NI 


了 调制 时 间 中 。 但 是 ， 旋 转 部 伯 
量 和 体积 ; 而 如 果 波 片 采用 切入 切 出 的 方式 ， 其 调 


F 的 使 用 会 导致 仪 费 拌 动 ， 从 而 降低 测量 精度 ， 增 加 系统 的 质 


出 时 间 将 大 大 增加 ,一 次 切换 时 间 一 般 大 


于 5 Fb. 这 些 缺 点 对 于 星 载 仪 器 更 为 严重 , 并 且 有 机 械 故 障 的 风险 。 相 对 而 言 , 波 片 +DKDP 


晶体 的 调制 方式 缩短 了 一 些 测量 时 间 , 但 
晶体 快速 调制 的 优点 ,仍然 需要 大 量 调 种 


日 于波 片 结构 的 存在 , 没有 能 够 充分 地 利用 DKDP 
时间, 实时 性 并 不 好 ,DKDP 晶体 和 液晶 材料 LCVR 


都 是 通过 改变 偏振 相位 的 方式 实现 对 入 射 光 的 偏振 调制 , 二 者 在 材料 、 性 能 等 方面 存在 较 大 


区 别 。 相 比较 而 言 ，LCVR 调制 频率 100Hz 左右 ， 对 温度 敏感 ， 光 谱 范 围 罕 ， 稳 定性 相对 


较 差 ， 在 实际 应 用 中 尚未 成 熟 。DKDP 材料 调制 频率 较 高 ， 可 达 1000Hz LA 


且 在 多 套 系统 中 成 功 应 用 ， 技 术 成 熟 可 靠 。 


基于 以 上 背景 ， 本 文 提出 


H DKDP 晶体 作为 偏振 调制 器 从 
而 DKDP 晶体 的 电压 -相位 延迟 响应 速度 可 以 达到 毫秒 级 


J ÆJ 


FEF， 重复 性 好 ， 


XX DKDP 晶体 的 实时 偏振 测量 方案 。 该 方法 采用 了 两 


8, RU 


,通过 外 加 电压 的 方式 来 控制 DKDP 晶体 的 相位 延迟 变化 ， 
加 方法 没有 任何 的 机 械 


运动 ， 没有 移动 部 件 ， 可 以 提供 写 秒 级 的 快速 切换 时 间 ， 再 加 上 曝光 时 间 短 ， 时 间 分 辨 紊 可 


达到 传统 波 片 +DKDP 调制 方式 的 10 倍 以 上 , 使 得 冻结 大 气 消 流 成 为 可 能 。 为 了 进 


步 验证 


该 方案 的 可 行 性 ， 本 文 在 方法 原理 、 参 数 优化 、 仿 真 分 析 和 测量 时 间 等 方面 进行 了 分 析 与 验 


证 。 


1 原理 和 方法 


光 的 偏振 信息 可 以 用 电 矢 量 、 琼 


斯 矢量 和 斯 托 克 斯 矢量 表示 。 其 中 电 矢 量 和 琢 斯 矢量 只 


能 对 完全 偏振 光 进 行 表示 ， 斯 托 克 斯 〈Stokes) 矢量 用 四 个 光 强 参量 可 以 实现 对 任意 偏振 态 


的 光 进 行 表示 ， 因 而 本 文选 用 


JFET 


体 表 达 式 如 下 : 


S = [L Q, U,V]? 


| =< E,E,' + EyEy* > Ey’ +E,’ 


Q-«E,E,' — E,E,' = E,?-Ey’#(1) 
U =< EEy + EyE,” >= 2E,Eycosó 


斯 矢量 来 表示 被 测 光 的 偏振 态 。 斯 托 克 斯 矢量 


S INR 
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V = -i< EE," + EJE,' >= 2E,Eysinó 


Jupes, 6 = yy — qu den X 与 Y 分 量 的 相位 差 : I 表示 光波 的 总 光 强 ，Q 表示 
将 光 分 解 到 X，Y 方向 后 ，X 与 Y 分 量 光 强 的 差 值 ，U 表示 将 光 分 解 到 45” 和 135° 方向 


Ja, 45° Fl 135° 分 量 光 强 的 差 值 ，Q，U 都 表示 光波 的 线 偏振 特 怕 


lim 


E. V 表示 将 光 分 解 为 右 


旋 和 左旋 分 量 后 ， 右 旋 与 左旋 分 量 光 强 的 产值 ，V 代表 了 光波 的 贺 
1.1 传统 波 片 与 DKDP 晶体 的 调制 方法 
传统 的 偏振 测量 方案 采用 了 4/4 波 片 与 DKDP 晶体 作为 偏振 调和 


偏振 片 透 光 轴 角 度 为 0” 。 


局 光 特 性 。 


Hart, 并且 一 次 完整 的 


斯 托 克 斯 矢量 测量 需要 进行 6 次 调制 得 到 6 组 光 强 值 . 此 方案 中 DKDP 蝇 体 快 轴 和 角度 为 45”， 


QWP DKDP Polarizer Detector 


图 1 传统 4/4 波 片 +DKDP 偏振 分 析 器 示意 图 


Figure 1 Schematic diagram of traditional 4/4 wave plate + DKDP polarization analyzer 


在 测量 纵向 磁场 时 ， 四 分 之 一 波 片 在 光路 之 外 ， 在 DKDP 蝇 体 上 加 和 矩形 正 负 高 压 ， 使 


得 相位 延迟 在 4/4 之 间 调 制 ， 左 右 圆 偏振 光 交 替 通 过 分 析 器 ，Stokes 矢量 V 信和 号 可 以 从 右 


旋 圆 偏振 光 与 左 旋 圆 偏振 光 的 差 中 得 到 。 测量 横向 磁场 时 ,将 四 分 之 一 波 片 放 在 光路 中 ,将 
其 光 轴 方向 角 调 为 0” ， 然 后 再 DKDP 晶体 外 加 矩形 正 负 高 压 ， 则 偏振 方位 角 分 别 为 0” 和 
90° 的 线 偏 振 光 交替 通过 检 偏 器 ， 它 们 的 差 值 即 Stokes 矢量 U 信号 。 调 整 四 分 之 一 波 片 的 
光 轴 方位 角 为 44”， 在 DKDP 昂 体外 加 矩形 正 负 高 压 ， 偏 振 方位 角 分 别 为 45” 和 135” 的 


线 偏振 光 交 蕉 通过 ， 相 减 后 得 到 Q fiim. de) res E] Vs hill BC 
表 1 4/4 波 片 +DKDP 偏振 调制 参数 


下 表 。 


Tablel 4/4wave plate +DKDP polarization modulation parameters 


A/4wave plate DKDP crystal Polarizer Measurement result 
45° -90° 0° I, =I+Q0 
45° 90° 0° I,=I-Q 
0° -90 0° Iz =I1+U 
0° 90° 0° I,=I-U 

90° 0° I, =I1+V 
-90° 0° 16 = 了 一 了 


在 此 方案 中 ， 四 分 之 一 波 片 的 切换 时 间 为 秒 量 级 ， 且 在 测量 


个 波长 点 时 需要 采集 6 


帧 光 强 图 像 ， 切 换 3 次 四 分 之 一 波 片 。 在 太阳 磁场 测量 中 ， 相 比 单 波长 点 测量 来 说 ， 多 波长 


点 测量 可 以 避免 受到 饱和 效应 造成 的 非 线性 、Doppler 位 移 以 及 Faraday 效应 等 的 影响 。 多 
波长 点 测量 一 般 选 取 谱 线 轮廓 上 的 6 个 波长 点 ， 此 时 ， 一 次 完整 的 谱 线 测量 至 少 需要 采集 
36 帧 光 强 图 像 ， 切 换 18 次 四 分 之 一 波 片 ， 导 致 测量 时 间 过 长 ， 时 间 分 辩 紊 过 低 ， 实 时 性 较 
Lm 
1.2 基于 双 DKDP 晶体 的 调制 方法 

当 光 的 偏振 态 用 斯 托 克 斯 矢量 S 表示 时 ， 传 输 介 质 对 光 偶 振 态 的 作用 可 以 用 称 勒 矩阵 
M 表示 。 光 束 依次 经 过 1…n 个 传输 介质 后 ， 斯 托 克 斯 矢量 为 : 

S’ = M4Mqa 477: M;M,SH(Q2) 


= 


偏振 测量 的 原理 如 下 图 : 


la 


Ma 
Polarization Mb Intensity lb 


modulator M, detector |, 
Mg la 


图 2 偏振 测量 原理 图 


Figure 2 Polarization measurement principle diagram 


偏振 测量 仪 通 常 由 偏振 调制 费 和 光 强 探测 器 构成 。 对 于 斯 托 克 斯 矢量 S (未 知 〉 的 待 测 


偏振 光 , 经 过 调制 状态 a 后 , 其 斯 托 元 斯 矢量 为 $4 = [la Qa Ua Va] ,此 时 调制 状态 a 对 应 的 
BIEREN Mao 


Q411 012 i3 Q14 


. 1921 422 423 24 #(3) 
Q31 032 433 34 


Q41 Q42 Q43 Q44 
Sq = M,S 
光 强 探测 器 只 能 探测 到 光 强 信息 ， 即 工 的 值 。 满 足 : 
la 2a41ILTa4,,Q + a13U + a,,V#(4) 
经 过 a，b，c，d 四 组 不 同 的 调制 状态 后 可 以 得 到 四 组 光 强 值 ， 满 足 ; 


Iq Q41 0412 Qiz Aya] [I 
Ip} |bu bui ps ao 
l| ca C2 C3 Ca FU 8G) 


la dij di diz dullv 
将 上 式 记 为 1% =XS， 当 调制 矩阵 X AY, TAME A, MUTE EI : 
S = X“11#(6) 
我 们 提出 了 一 种 基于 双 DKDP 的 偏振 测量 方法 ， 如 图 3 所 示 。 该 偏振 分 析 器 的 偏振 调 
制 器 由 两 组 快 轴 方 向 不 同 的 DKDP 晶体 和 一 个 偏振 片 构成 ， 光 强 探测 器 为 光 强 采集 相机 。 
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DKDP1  DKDP2 Polarizer Detector 


图 3 基于 DKDP 的 偏振 分 析 器 结构 示意 图 


Figure 3 Schematic diagram of the structure of the polarization analyzer based on DKDP 


假定 入 射 光 斯 托 克 斯 矢量 为 Sim， 


I 
Q 
U 
V 
则 通过 偏振 调制 器 后 的 其 斯 托 克 斯 矢量 变 为 Su， 

Sout = My - M2 : M1 - Sin#(8) 


#(7) 


Sin = 


其 中 ，Mi，M2，Mp 分 别 为 第 一 组 DKDP 品 体 、 第 二 组 DKDP di A RU f as A PAE 


4 DKDPI 晶体 快 轴 角 度 为 0” ， 相 位 延迟 量 为 8，DKDP2 晶体 快 轴 角 度 为 43” ， 相 位 
延迟 量 为 g， 检 偏 器 〈 偏 振 片 ) 透 光 轴 为 0” 。 检 偏 器 的 偏振 方向 角 与 笛 卡 尔 坐 标 系 的 选取 
有 关 ， 为 了 方便 分 析 ， 我 们 此 处 令 检 偏 器 的 偏振 方向 为 坐标 系 的 X 轴 方 向 。 则 DKDP 晶体 


和 检 偏 器 的 穆 勒 矩阵 Mi ，M2，Mo 为 : 


1 0 0 0 
..[0 1 0 0 
ur 0 0 cosó a #(9) 
0 0 -—sinó cosó 
1 0 0 0 
M. = 0 cose 0 #(10) 
2~ 10 0 1 0 
0 sine 0 cose 
1100 
M, = 0.5 | A :| #(11) 
0.00 0 
此 时 ， 通 过 偏振 分 析 器 后 的 斯 托 克 斯 矢量 为 ; 
I 
Sout = = Mp M5: My : Sin# (12) 
V 


光 强 探测 器 接收 到 的 光 强 为 : 
I = 0.5( + Qcose + Usinesin6 — Vsinecos6)#(13) 
改变 两 组 DKDP 晶体 的 外 加 电压 ， 得 到 4 组 不 同 的 偏振 调制 状态 ， 光 强 探测 器 测 得 4 
组 不 同 的 光 强 值 10wt = [Lilo Is, I4]T 。 由 于 一 般 只 需要 斯 托 克 斯 矢量 的 相对 值 ， 所 以 式 中 的 
0.5 可 略 去 ， 得 到 偏振 调制 矩阵 为 : 
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cose,  sing4sinóÓ, —sine,cosd, 
COSE; sinez,sind, —sine,cosd, 
COSE, sinessind;  —sin£4CosÓ, 
COSE, sine,sind, —sine,cosd, 


WAT out = 0.5X .Su ， 当 偏振 调制 矩阵 X n3, nf 


BRR 


E 


I I 
_ [Q| _ >ya |2 
Sin = = 2X #(15 
mn U Ls ( ) 
V n 


Sin 即 为 待 测 光 斯 托 克 斯 矢量 。 


2 参数 最 优化 设计 


vT 偏振 测量 过 程 中 测量 精度 会 受到 偏振 调制 器 调制 矩阵 误差 和 光 强 探测 器 探测 误差 等 的 
人 影响 ,而 偏振 调制 矩阵 参数 的 不 同 选 择 对 于 误差 的 传播 影响 不 同 。 在 保证 偏振 调制 矩阵 可 逆 
e 的 同时 , 合理 地 选用 调制 参数 可 以 显著 减 小 噪声 扰动 和 系统 误差 对 测量 结果 的 影响 。 此 处 以 
假定 只 存在 光 强 探测 器 误差 为 例 。 
AV 为 光 强 探测 值 ，I 为 光 强 真 值 ，S 为 斯 托 殉 斯 矢量 计算 值 ，S 为 真 值 。 光 强 误差 A1 为 
AI = 了 一 17， 斯 托 克 斯 矢量 误差 AS 为 
AS-S'—S-XI'—S-—X^(IT AD) - Se XA 


AI 由 1 = XS， 得 
; IHIXSI < IXS 
I|ASI| = X74 Ar || < XHAI) 


IASI y—1 lA 
LX xy 
ise MIX T 


上 式 给 出 了 相对 误差 |IAsIMllsll 的 误差 上 界 ， 光 强 探 测 器 误差 在 求解 过 程 中 可 能 会 被 放大 
1 一 INIXL 倍 ， 俩 振 调 制 矩 阵 误 差 结 果 相 同 ， 在 此 不 做 推导 。 通 常 称 
cond(X)-l X HIX AERE X 的 条 件数 CCN) PP. AE E X 的 条 件数 越 小 ， 在 矩阵 方程 求解 过 
程 中 误差 的 放大 越 小 。 根据 研 究 结 果 以 及 仿真 验证 ,对 于 全 斯 托 克 斯 矢量 测量 , 偏振 调制 甜 
阵 的 条 件数 最 小 值 为 V3， 因 此 ， 当 调制 矩阵 条 件数 为 V3 时 ， 偏 振 调制 参数 最 优 。 
偏振 调制 矩阵 如 式 〈14) 所 示 ， 基 于 优化 速度 更 快 的 拟 牛 顿 算法 ， 通 过 最 小 化 矩阵 X 
的 CN 来 实现 优化 过 程 5"。 拟 牛顿 算法 是 解决 无 约束 优化 问题 有 效 且 成 熟 的 算法 ,相对 于 牛 
顿 算法 ， 其 只 要 求 目标 函数 一 次 连续 可 微 ， 易 于 实现 ， 计 算 量 小 。 通 常情 况 下 ， 实 现 相 同 的 
最 佳 优化 CN 有 无 限 种 解决 方案 ， 用 于 偏振 测量 时 ,它们 具有 相同 的 性 能 ， 此 处 我 们 仅 选取 
其 中 一 组 参数 : 


表 2 偏振 调制 方案 参数 


Table 2 parameters of polarization modulation scheme 


DKDPI delay (rad) DKDP2 delay (rad) 


Polarizer angle 


5.4978 5.3279 0° 
5.4978 0.9553 0° 
3.9270 2.1863 0° 
3.9270 4.0969 0° 
此 时 ， 偏 振 调制 矩阵 为 
1 0.5773 0.5774 0.5774 
y,z-2|i 05773 —0.5774 —0.5774 
cN |1 —0.5773 —0.5774 0.5774 
1 —0.5773 0.5774  —0.5774 
此 时 ， 偏 振 调制 矩阵 求解 时 对 误差 的 放大 最 小 。 


Wt» 


， 偏 振 调制 矩阵 的 等 加 权 方 差 (EWV) 也 是 评估 误差 传递 的 标准 之 一 ， 但 EWV 


的 值 主要 受到 数据 见 余 的 影响 。 并 


， 本 文中 的 新 测量 方案 仅 适 ) 


] 于 带宽 极 罕 (AX<0.1A) 


应 用 场景 下 ,通常 与 滤 光 器 配合 使 用 ， 其 光谱 色散 问题 可 以 忽略 , 不 考虑 偏振 测量 效率 。 本 


文 所 使 用 的 DKDP 晶体 透 过 范 


3 仿真 分 析 验 证 


在 传统 的 波 片 +DKDP 蝇 体 式 偏振 分 析 器 中 ， 常 采 
的 存在 ， 其 精度 会 略 优 于 同等 条 件 下 CCN=V3) 的 4 次 调制 方案 。 


A 


围 为 200~2150nm， 在 目标 波长 532.4nm 处 透 过 率 可 达 9796. 


j 6 次 调制 的 测量 方案 ， 由 于 宛 余 项 


在 此 我 们 比较 两 者 的 精 


度 差异 ， 并 验证 提出 的 DKDP 侦 振 分 析 器 精度 是 否 符合 测量 要 求 。 


在 本 文中 太阳 偏振 测量 


精度 要 求 为 2Xx10-3。 由 于 市 面 上 常见 优 
误差 可 以 忽 
及 DKDP 晶体 和 偏振 片 角度 误差 对 测量 结果 影响 。 

1) 光 强 探测 器 误差 


测量 过 程 中 , 受到 环境 杂 散 光 和 光 强 探测 器 本 喘 噪声 的 影响 , 光 强 探测 器 测 得 的 光 强 人 


ij 振 片 消光 比 误差 约 为 10- 一 10， 其 引起 的 测量 
各 ; 在 此 ， 我 们 将 比较 两 种 调制 方案 中 光 强 探测 器 误差 、 


DKDP 相位 延迟 误差 以 


存在 一 定 的 误差 。 此 误差 将 直接 反应 到 竺 测 光 的 偏振 误差 ， 由 此 我 们 有 必要 验证 提出 测量 方 


案 的 测量 误差 满足 测量 精度 。 


四 


= pH SP 


一 般 情 况 下 ， 光 强 探测 器 的 非 线 度 在 1-296, XEULi AL fi Hel fg 


响应 函数 进行 标定 。 假 定 光 强 探测 器 的 非 线 性 误差 为 5=0.02， 
强 误差 导致 的 偏振 误差 计算 方案 如 下 
AS = S'—S = X7]' — X-11#(17) 


为 了 评估 方案 的 性 能 ， 
入 射 偏振 态 ( 不 同 偏振 分 量 ) 的 测量 精度 。 


ea 


(Poincaré Sphere) 球面 上 的 每 一 个 点 来 表示 。 我 们 选取 Poincaré BRI 


点 ， 其 中 20 个 不 同 椭圆 


Hz 


通常 都 需要 对 其 


且 只 考虑 光 强 探测 器 误差 , 光 


我 们 分 析 了 光 强 噪声 为 0.02 和 0.001〈 标 准 差 ) 时 ，1000 个 不 同 
ABH, TER HI TEE Ah He AS AT DA 


J PE DN SE BR 


fl Ef 1000 个 均匀 采样 


度 的 的 圆 从 南极 到 北极 ,每 个 圆 的 50 个 方位 角 顺 时 针 方向 移动 指向 


的 偏振 态 作为 入 射 偏振 态 。 在 仿真 中 ， 通 过 添加 一 个 信号 独立 的 高 斯 噪声 ， 每 个 状态 产生 
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100 组 强度 值 ， 计 算出 100 组 测量 误差 的 标准 差 作 为 对 应 偏振 态 的 测量 误差 。 


显然 ， 两 种 测 


m. 
0.5 8 = 
D 
e 0 T x 
[2] n 
[^ 
as 
-0.5 B 
N 
4 i. 
1 EE 
0.5 1 S. 
0 0.5 es 
-0.5 05 0 REN" m —»» oOo mw o m 9 
s2 a 4 $1 Polarization State 
图 4 庞 加 莱 球 图 5 偏振 采样 点 
Figure 4 Poincare sphere Figure 5 polarization sampling points 


四 


6 中 , Q1-V1 为 我 们 提出 的 DKDP 偏振 测量 方案 结果 , Q2-V2 为 传统 方案 测量 结果 。 


EL 


量 方案 的 斯 托 殉 斯 天 量 误差 在 同一 数量 级 ， 但 由 于 光 强 探测 器 非 线性 在 2%, 


导致 了 两 种 方案 的 误差 都 达到 了 10- 误 差 等 级 ， 难 以 满足 测量 要 求 。 我 们 又 将 光 强 探测 器 
的 非 线性 降低 到 1%o 后 ， 得 出 了 图 7 中 结果 。 可 以 看 到 ， 在 光 强 探测 误差 为 1% 时 ， 对 于 斯 
托 克 斯 拓 量 Q U. V 分 量 的 测量 误差 ,传统 测量 方案 的 误差 在 1.4 x 10 飞 左右， 而 新 的 偏 
振 测量 方案 误差 在 1.7 x 10 飞 左右， 误差 均 小 于 2 x 10 飞 ， 满 足 偏振 测量 要 求 。 


Standard Deviation 


x103 
r 


Ee 


Standard Deviation 
5 


i 


i 10 p 300 p 500 EE 700 ET 900 1000 ý 0 100 E ao 400 300 600 700 300 ms 1000 
Polarization state (670.02) Polarization state (670.001) 
图 6 Stokes 矢量 误差 比较 (87296) 图 7 Stokes 矢量 误差 比较 (6=1%o) 
Figure 6 Stokes vector error comparison (6=2%) Figure 7 Stokes vector error comparison (6 = 1 %o) 
2) DKDP 品 体 相位 延迟 误差 


测量 过 程 中 ， 在 控制 外 加 电压 切换 DKDP 品 体 相位 延迟 时 ， 由 于 电压 -相位 延迟 曲线 会 


存在 一 定 的 误差 ， 导致 偏振 测量 的 结果 也 会 出 现 误差 。 经 过 对 应 波长 下 (本 文 为 532.4nm ) 


的 电压 -相位 延迟 标定 后 ，DKDP 唱 体 的 延迟 误差 一 般 小 于 0.001 弧度 。 我 们 需要 验证 两 种 


测 | 


a TP H 


日 相位 延迟 误差 引起 的 偏振 测量 误差 是 否 在 误差 允许 范围 内 。 被 测 偏 振 态 的 采 


样 方法 与 上 述 方法 相同 。 在 仿真 时 ， 添 加 一 个 独立 的 标准 差 为 0.001 的 角度 高 斯 噪声 ， 得 到 
1000 个 偏振 状态 的 测量 误差 。 
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Standard Deviation 


0 1 1 1 1 1 1 1 1 1 
0 100 200 300 400 500 600 700 800 900 1000 
Polarization state 


图 8 DKDP 相位 延迟 误差 导致 的 测量 误差 比较 


Figure 8 Comparison of measurement errors caused by DKDP phase delay error 

8 P, QI-VI 为 我 们 提出 的 双 DKDP 偏振 测量 方案 结果 ，Q2-V2 为 传统 方案 测量 结 
果 。 容 易 看 出 ， 两 种 测量 方式 由 DKDP 相位 延迟 误差 引起 的 偏振 测量 误差 均 小 于 10-3， 其 
中 QQ 分 量 测 量 结果 优 于 传统 测量 方式 ，U 分 量 和 V. 分 量 误 差 值 几乎 相等 ， 满 足 测量 需求 。 

3) DKDP 晶体 快 轴 与 偏振 片 角度 误差 

由 于 DKDP 晶体 与 偏振 片 在 偏振 分 析 器 中 国定 不 动 ， 在 对 DKDP 晶体 的 快 轴 与 偏振 片 
角度 完成 标定 后 , 其 角度 误差 远 小 于 0.001 弧度 。 此 处 我 们 假定 DKDP 晶体 的 快 轴 与 偏振 片 
角度 误差 均 为 0.001 弧度 , 验证 其 是 否 满足 测量 精度 的 要 求 。 使 用 相同 的 采样 方法 得 到 1000 
个 偏振 采样 点 ， 并 添加 一 个 独立 的 标准 差 为 0.001 的 高 斯 噪声 到 DKDP2 晶体 的 快 轴 角度 参 
Be, 得 到 受 快 轴 角 度 误差 影响 的 偏振 测量 误差 。 同 理 得 到 偏振 片 角 度 误差 导致 的 偏振 测量 误 
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图 9 DKDP2 快 轴 角度 误差 结果 图 10 偏振 片 角度 误差 结果 
Figure 9 DKDP 2 fast axis angle error results Figure 10 polarizer angle error results 


从 图 9 和 图 10 中 我 们 可 以 看 出 ， 在 两 种 角度 误差 的 影响 下 ， 测 得 的 斯 托 克 斯 矢量 中 只 
有 Q 分 量 会 受到 DKDP2 晶体 快 轴 角 度 误 差 和 偏振 片 角度 误差 的 影响 ， 且 其 测量 误差 小 于 
2 x 10-3， 满 足 偏振 测量 的 要 求 。 由 于 DKDP 晶体 固定 不 动 ， 受 到 快 轴 角 度 误 差 的 影响 远 远 
小 于 传统 的 波 片 式 结构 。 
最 后 ,我 们 对 两 种 方案 的 测量 时 间 进 行 了 对 比分 析 。 其 中 , DKDP 晶体 毫秒 级 的 响应 时 


间 可 忽略 不 计 ， 令 相机 采集 频率 为 10 帧 / 秒 ， 为 了 提高 精度 ， 每 个 调制 位 置 采集 两 帧 图 


表 3 方案 对 比 〈 单 波长 点 ) 


Table 3 scheme comparison (single wavelength point) 


像 。 


Modulation Modulation Camera image frames Measuring time 
frequency steps frequency 
Traditional scheme >5s 3 (wave plate) 10 frames / 12 frames >16.2s 
sec 
Double DKDP millisecond 4 10 frames / 8 frames ~0.8s 
scheme sec 
364 方案 对 比 〈6 个 波长 点 ) 
Table 3 Scheme comparison (6 wavelength points) 
Modulation Modulation Wavelength point image Measuring 
frequency steps switching time frames time 
Traditional >s 18 14s 72 >111.2s 
scheme 
Double DKDP millisecond 24 14s 48 ~18.8s 
scheme 
从 上 面 的 对 比 中 我 们 可 以 看 出 ， 在 单 波长 点 偏振 测量 中 ， 传 统 方案 测量 时 间 大 于 16.2s 


E EM 


EXN DKDP 测量 方案 测量 时 间 〈0.8s) 的 20 fis 


以 上 ; dE 6 个 波长 点 


的 偏振 测量 中 ， 传 统 方 


HB. 


时 间 约 为 18.8s， 由 
影响 存在 ， 传 统 方案 测量 时 间 为 双 DKDP 方案 的 6 倍 左右 。 


案 测 量 时 间 大 于 111.2s, X DKDP 测量 方案 测 
的 
XX DKDP 晶体 原理 的 实时 1 
可 分 别提 升 20 倍 和 6 fit. 


里 


= 


于 波长 点 切换 频率 较 低 
因此 ， 采 用 本 文 提出 的 基于 
有 振 测 量 方案 ， 在 单 波长 点 和 多 波长 点 偏振 测量 中 ， 时 间 分 辨 率 


太阳 大 气 探测 在 空间 天 气 预报 和 太阳 物理 研究 等 方 


生发 挥 着 重要 作用 , 而 太阳 磁场 驱动 


太阳 活动 爆发 的 根本 驱动 力 , 偏振 测量 是 实现 太阳 人 磁场 探测 的 重要 过 程 。 目前 偏振 测量 的 主 


要 手段 包括 旋转 波 片 调制 方式 与 波 片 +DKDP 晶体 调制 方式 ， 
相 比 旋转 波 片 可 以 减少 旋转 波 片 的 次 数 ， 缩短 测 
大 量 时 间 ， 实 时 性 不 佳 ， 
速度 ， 本 文 提 出 了 基于 双 DKDP 


晶体 原理 的 实时 偏振 测量 方法 ， 采 用 


波 片 +DKDP 唱 体 的 测量 方法 
量 时 间 , 但 由 于 波 片 结构 的 存在 ,仍然 需要 
容易 出 现 仪 器 拌 动 等 问题 。 基 于 DKDP 晶体 毫秒 级 的 电光 响应 
了 两 组 DKDP 


H 
HH 


体 作 


为 俩 振 调制 器 件 ， 调 制 时 间 可 从 数秒 提升 到 毫秒 级 ， 大 大 提升 了 测量 的 实时 性 。 本 文通 过 仿 


真 验证 了 该 方法 的 可 行 性 ， 
间 分 辩 率 可 提升 全 传统 方案 的 20 倍 和 6 fi, OSEE 


在 精度 上 优 于 2 x 10-， 在 单 波长 点 和 多 波长 点 偏振 测量 中 ， 时 
大 气 实时 偏振 探测 有 


E- e N 
EE 要 意义 。 
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Abstract: Magnetic field is the fundamental driving force of solar eruptions. High precision 
real-time polarization measurement is a common method to realize magnetic field measurement. 
However, the commonly used wave plate measurement method needs to be measured many times, 
and there is a mechanical rotation structure, which is easy to cause instrument jitter and reduce the 
measurement accuracy; the rotation or switching of wave plate will also lead to a longer 
measurement time. At present, the measurement period of solar atmospheric polarization based on 
DKDP crystal and wave plate is long, and the real-time performance is not well, so it can not be 
applied to the detection of rapidly changing solar active magnetic field. Based on the above 
background, this paper proposes a real-time polarization detection technology of the solar 
atmosphere based on the principle of double DKDP crystals. Two DKDP crystals with different 
fast axis directions are used to realize the polarization modulation of the incident sunlight and 
obtain the polarization Stokes vector of the solar atmosphere. Compared with the traditional 
method of combining rotating wave plate and DKDP crystal, the proposed method can increase 
the single polarization modulation rate from several seconds to millisecond, which greatly 
improves the real-time performance of solar magnetic field measurement. The simulation results 
show that the measurement accuracy of the proposed method is 2 x 1073, and the measurement 
time resolution of single wavelength point is improved by more than 20 times, which can realize 

the real-time accurate measurement of the polarization information of solar magnetic field. 
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